Spin Injection in a Ballistic Two-Dimensional Electron Gas 
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We explore electrically injected, spin polarized transport in a ballistic two-dimensional electron 
gas. We augment the Biittiker-Landauer picture with a simple, but realistic model for spin-selective 
contacts to describe multimode reservoir-to-reservoir transport of ballistic spin 1/2 particles. Clear 
and unambiguous signatures of spin transport are established in this regime, for the simplest mea- 
surement configuration that demonstrates them directly. These new effects originate from spin 
precession of ballistic carriers; they exhibit strong dependence upon device geometry and vanish in 
the diffusive limit. Our results have important implications for prospective "spin transistor" devices. 

PACS: 71.70.Ej, 73.40.-c, 85.90.-h, 73.50.-h 
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The concept of a spin transistor, first proposed almost 
ten years ago has attracted widespread interest ||] 
but its experimental realization remains elusive 0. It 
is based upon electrical injection of spin polarized carri- 
ers from a ferromagnetic conductor into an electron gas 
within a semiconductor. Electrons propagating in the in- 
terfacial electric field confining them to the device chan- 
nel experience an effective magnetic field that induces 
spin precession; this is called the Rashba effect Q{. The 
rate of spin precession should be tunable through an ex- 
ternal gate voltage, which will add to the confinement 
potential |||(| With ferromagnetic source and drain con- 
tacts acting as spin polarizer and analyzer, an electron 
device analogous to an electro-optic modulator is en- 
visaged. 

Experiments to date demonstrate that the spin transis- 
tor geometry [Fig. |l|(a)], in general, leads to strong Hall 
phenomena that are unrelated to true spin transport || . 
They arise due to the requisite proximity of miniature 
magnets and the low density (and, hence, high Hall co- 
efficient) electron gas. Since these Hall phenomena de- 
pend directly upon the magnetization state of these mag- 
netic contacts, they often closely mimic the signals ex- 
pected from spin transport experiments — especially those 
in which the relative magnetic orientation of the "spin po- 
larizer" and "analyzer" contacts is varied. However, in 
early experiments on (diffusive) spin injection in met- 
als spin precession phenomena provided an alter- 
nate and crucial experimental proof of spin transport. 
In this Letter, we establish for the first time the anal- 
ogous, and unambiguous, experimental signatures to be 
expected from spin injection and precession phenomena 
in a ballistic two-dimensional electron gas (2DEG). The 
observation of these precessional effects will constitute a 
definitive experimental demonstration of electrical spin 
injection in semiconductor systems. 

The spin transresistance, Rs [Fig. 0(a)] provides the 
most direct demonstration of spin transport This 
nonlocal transport coefficient is free of obfuscating back- 
ground signals unrelated to spin injection. In the diffu- 
sive limit, if current contact Fl is replaced with one that 
is unpolarized, no voltage will appear between the ana- 



lyzer contact F2 and a suitably defined ground reference 
R. In this case these voltage contacts, being well out- 
side the net current path, remain at equipotential. With 
a polarized current contact Fl, injected magnetization 
can lead to steady-state spin accumulation that persists 
over the entire length of the channel if 5 s ^ L. This in- 
duces disequilibrium between the local spin-resolved elec- 
trochemical potenti als of F 2 and R, and yields a finite 
R s - Here 6 S = ^JIqIs]2 is the spin diffusion length, 
is = vft~s and £q = vftq are the spin and momentum 
mean free paths, t§ and tq arc the effective spin and mo- 
mentum relaxation times, and vf is the Fermi velocity. 

In the ballistic regime, however, it is not appropri- 
ate to speak of spin accumulation since a local chemi- 
cal potential cannot be meaningfully defined within the 
channel. Accordingly, our description of the ballistic 
spin transresistance is based upon Biittiker's picture for 
mesoscopic transport within a multiprobe conductor ||. 
Here we augment this with a new model describing spin- 
selective contacts. Our procedure is as follows: (a) 
We first develop a simple description of spin-selective 
contacts, based upon careful consideration of the fer- 
romagnetic/semiconductor (F/S) contacts in (our) real 
devices ||. (b) We construct an 8-reservoir model, af- 
ter Biittiker, to describe the spin injection experiment, 
(c) Boundary conditions are used to constrain the spin- 
resolved currents and chemical potentials. These lead 
to a simpler 4-reservoir problem for the spin transresis- 
tance, in terms of reservoir-to-reservoir, spin-resolved 
transmission probabilities, T^, of the 2DEG forming the 
device conduction channel. Here the indices i, j specify 
the reservoirs themselves, and a, f3 their constituent spin 
bands, (d) The requisite are then calculated semi- 
classically, using a modified Monte Carlo numerical tech- 
nique (described below). We follow the electrons' ballistic 
trajectories and the phase of their spin wavefunctions as 
they pass through the device, while ignoring the phase of 
their spatial wavefunctions. For unpolarized ballistic sys- 
tems, this semiclassical approach has proven remarkably 
consistent with experimental data at T ~ 4K, where the 
electron phase coherence length is smaller than typical 
dimensions of nanoscale devices 110(1 . 
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FIG. 1. Model for ballistic spin injection in 2D. (a) Mea- 
surement configuration: a current, /, is injected through the 
2DEG via a ferromagnetic contact Fl and an ohmic contact 
L. The spin transresistance, Rs = V/I, arises from spin po- 
larized carriers traversing a distance L from the path of the 
net current, which induce a nonlocal voltage, V, between a 
second, similar, pair of contacts, F2 and R. (b) The conduc- 
tor beneath the ferromagnetic contacts (DSPR) is assumed 
to be a Disordered, but Spin Preserving Region, (c) The full 
8-reservoir model; complete ellipses represent spin-relaxing 
reservoirs, half ellipses represent spin-resolved reservoirs. Fl 
and L are current contacts, F2 and R are voltage probes. 
T"? denotes the 2DEG device channel in which spin preces- 
sion occurs. Other multimode leads are denoted by three 
arrows and ellipsis. Panels (d), (e), (f) illustrate decomposi- 
tion of the 8-reservoir model. (Panel (e) depicts the reduced 
4-reservoir problem.) 



ing physics. Of course, many complexities in real devices 
may diminish spin transport effects fl2|| . Here our aim 
is to establish what may be expected in ballistic systems 
under optimal conditions. 

Measurement of Rs involves four terminals [Fig. 0(a)], 
two that are spin-selective, Fl and F2, and two that 
are conventional, i.e. momentum- and spin-relaxing, L 
and R. As depicted in Figures 0(d, e, f), the full prob- 
lem separates into three sub-components. Fig. 0(d) rep- 
resents the spin-up and spin-down currents (/if, hi) 
that flow between Fl, If, If, and L. A Sharvin resis- 
tance HI, R sh = (h/2e 2 )(k F w)/n = (h/2e 2 )N ch , arises 
between If, If and the multichannel conductors con- 
necting them to L. Under conditions of current flow this 
yields the spin-resolved electrochemical potential differ- 
ences /iif - n L = 2eR sh lL] and /z u - fi L = ZeRshhi- 
Here, the factors of 2 arise because transport is spin 
resolved; hp, w, and N c h are the Fermi wave vector, 
channel width, and number of occupied modes within 
the 2DEG device channel, respectively. Similarly, at 
the rightmost side of Fig. 0(f), current flow between 
the reservoirs 2f , 2f and R establishes the electrochem- 
ical potential differences fi2i — t l R — ^&Rshh] and 
M21 - M-R = ZeRshlin- Also, [i F 2] = M2T smce 110 cur- 
rent flows between these reservoirs. Note that all Ps 
here represent net currents (forward minus reverse con- 
tributions). In our model, the following sum rules hold: 
I = hi + hh I = II] + hi, = Jr| + I Rl , and 
^lf + hi = h\ + hi = 0. As the reservoirs in Fig. 0(f) 
are voltage contacts, net current is conserved separately 
for each spin band, 1^ + h] = hi + hi = 0. These 
expressions can be manipulated to yield 



Figure 0(b) depicts our model for the spin-selective 
contacts which comprises two elements: F2, a fully 
spin-polarized reservoir which is in perfect contact 
with DSPR, a second Disordered (i.e. momentum- 
randomizing) but Spin- Preserving Region that consists 
of separate spin-up and spin-down bands. The separate 
spin-resolved reservoirs comprising the DSPR's (If, If, 
2f, 2f) model low mobility regions always present be- 
neath unalloyed ferromagnetic metal contacts in typical 
InAs devices ||. Disorder within them yields significant 
momentum randomization and, hence, a short £q. How- 
ever, in contrast to the usual picture describing unpo- 
larized reservoirs ||, we assume these special contacts 
are small compared to 5s, thus any spin disequilibrium 
within them is preserved. This, in fact, is consistent with 
the more restrictive constraint, 6s <] L, which is generic 
and fundamental to any spin injection experiment. If 
significant spin relaxation occurs anywhere in the device, 
including the vicinity of the ferromagnetic contacts, spin- 
selective transport is suppressed O]. In this Letter, for 
sake of clarity, we consider the most ideal situation, ini- 
tially assuming that Fl and F2 are fully polarized at the 
Fermi surface (half- metals) . This approximation serves 
to illustrate the most important aspects of the underly- 
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Given these relations, calculation of Rs reduces to a 
4-terminal problem that solely involves the four spin- 
resolved reservoirs: If, If, 2f, and 2f and the 2DEG 
device channel that connects them [Fig. 0(e)]. Modifying 
Biittiker's formula to account for the spin-resolved chan- 
nels, the 4-terminal linear response at zero temperature 
becomes 



(N ch - R%*)»ia 



rriQ.Q 



( 2 ) 



Transport within the ballistic multimode 2DEG con- 
ductor is fully represented by the transmission and reflec- 
tion coefficients, ly and Rf". These describe carriers 
incident from the lead i with spin polarization a, that are 
transmitted into lead j with final spin state (3 ; and car- 
riers incident from ia that are reflected back into same 
lead and spin channel, respectively. The coefficients f/y 

in Eq. satisfy the sum rule £\ Q ^iL. = U?f = 
ensuring the current sum rules of Eq.jll and that all cur- 
rents vanish when the \Xi are equal. 
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FIG. 2. Ballistic spin transresistance in an external field 
normalized to Bq = pp/ew, at which the cyclotron radius 
equals the channel width, (a) For a channel with L/w — 15 in 
perpendicular field, we plot two traces representing ujp /u) c = 1 
and 0.19, appropriate for a typical metal and for InAs, re- 
spectively. (b,c) Spin transresistance for three different con- 
figurations and two channel lengths, L/w = 3 and 15. Here 
u l /uj c = 0.19 (InAs). 



Simplification of Eq. [T] and 
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where p, L>R = n L . R /2eR sh and S = (1 + Uy l U . The 
elements of S satisfy the same sum rules that constrain 
U (for identical reasons). For parallel alignment of po- 
larizer and analyzer, Fl and F2, which we denote by the 
superscript (11), these steps yield 
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For antiparallel alignment, only the sign changes: 

We obtain the requisite elements of S numerically, 
extending the semiclassical billiard model |[o| to al- 
low tracking of an electron's spin wavefunction along 
ballistic trajectories linking the spin-resolved reservoirs 
(lt,l|,2t,2J,) at either end of the 2DEG device channel. 
We consider electrons confined within a hard- wall chan- 
nel, of length L and width w. The are calculated by 
injecting and following a large number of electron trajec- 
tories (typically > 10 4 ) propagating at vp |l6|| . 



For each path segment traversed by the electron be- 
tween boundary reflections, the phase of its spin wave- 
function evolves continously via the local Larmor fre- 
quency ujp — g*eB/2m. Here g* is the effective elec- 
tron g-factor, B the local magnetic field, and to the free 
electron mass. Total precession is accumulated for each 
complete trajectory, which is the sum of these segments. 
For each segment the electron's spin precession is calcu- 
lated analytically [jl4j and incorporated into the Monte 
Carlo procedure. 

In Fig. |5|(a) we display Rg as a function of per- 
pendicular magnetic field strength. The prominent and 
striking new feature is that R$ is oscillatory, a ballis- 
tic phenomenon not found in the diffusive regime. In 
Fig. ||(b,c) we display Rg calculated for three orien- 
tations of the external field — two that are in-plane and 
the perpendicular case, displayed again for comparison. 
In all three cases the Fl, F2 magnetizations are parallel 
and y-oriented. 

When the external field is along y, the injected carriers 
remain in spin eigenstates and do not precess. In this sit- 
uation Rg^ is a positive constant [Fig. ||(b,c)]. However, 
with an i-oriented field precession is maximal, and Rs 
oscillates. Since orbital effects are absent for an in-plane 
field, the oscillations in this case arise purely from spin 
precession and the oscillation period, AB, is determined 
by the condition 2-7m = oj^tr, *- e -, AB = h/ (3*MstTi?)- 
Here t^R — S/vp is a typical transit time from 1 — * 2, 
and [Lb is the electronic Bohr magneton. AB is thus 
inversely proportional to S, a typical path length aver- 
aged over the injection distribution function. The decay 
of Rs occurs on a field scale where ajpStpR ~ 7r; i.e. for 
B = hir I \g* lib StpR) beyond which precession amongst 
the different contributing trajectories tends to get out of 
step. Here StpR. = [tj-R ~ {^tr) 2 ] 1 ^ 2 is the variance in 
path lengths traversed while propagating from 1 — > 2. 

Perpendicular field (B ext \\z) is special — it induces both 
spin and orbital effects. (The characteristic field scale 
for the latter is Bq = pp/ew, at which the cyclotron ra- 
dius, r c ~ vf/lu c , equals the channel width, w.) The 
frequency ratio, lul/lu c — (g* /2)(m* /m) describes the 
relative importance of orbital and spin transport phe- 
nomena. Here, pp is the Fermi momentum, lo c — eB/m* 
the cyclotron frequency, and m* the effective mass. For 
InAs (to* = 0.025, g* = 15) this ratio is - 0.19, for In- 
GaAs ~ 0.1, whereas it is roughly 1.0 for most metals. In 
the latter spin and orbital effects have similar periodicity 
so disentangling them is difficult [Fig. ||(a)]. 

As mentioned, electrons confined within an InAs het- 
erostructure are subject to an internal Rashba field, 
present even for zero applied magnetic field. This can 
be modeled by a Hamiltonian [[|, Hp = a so [er x k] • z. 
Comparing Hp to the Zeeman term we write the ef- 
fective Rashba field as — 2a so k x z/(g*LiB). Here 
a so = Ap/2kp is the spin-orbit coupling parameter (A# 
is the Rashba splitting Jl5|]), and k and kp are the elec- 
tron and Fermi wave vectors, respectively. Using data 
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FIG. 3. Spin transresistance vs. reduced Rashba frequency, 
lur = 2m*a BO w/h 2 , at zero applied (external) field, for two 
different device channel lengths, L/w — 3 and 15. The pa- 
rameter u)r can be controlled by an external gate voltage. 
Shaded regions delineate the range of tunability expected for 
InGaAs devices |g| of three widths, 0.1 , 0.3, and 0.5 /im. 

from Heida et al. |ll| , we estimate this internal field to be 
about 5T for an InAs 2DEG. Since B/j is always in-plane, 
all electron trajectories are straight when the external 
field has no out-of-plane component. For fully polarized 
injection this yields the simple expression 

(TT) 1 — 2f 

Rs = (l + 2t)(2t-3) Rah ' (5) 

where t represents T^^i), normalized by N c h- 

Figure |^ displays how Rashba-induced spin precession 
is manifested in Rs for zero external field. We represent 
the effective Rashba field strength by the dimensionless 
frequency ujr = 2a so m*w/h 2 ; at ujr = 1 an electron pro- 
cesses one radian after traversing a distance w. As shown, 
the oscillations decay quickly initially, but exceedingly 
slowly thereafter. No spin precession occurs for lur = 0, 
hence, t = 1 yielding R^^ = B§^*~^ — R s h/3, a sim- 
ple result of current division. For finite wjj, Rs displays 
strong dependence upon the orientation of the magneti- 
zations M (of Fl, F2; assumed parallel), in relation to 
the device channel's principal axis (£). For M||x (par- 
allel to the channel), precessional effects are maximal. 
With increasing Rashba field, the variance in contribut- 
ing path lengths causes the oscillations in Rs to decay, as 
described previously for the case of finite external field. 
Here, however, the contributions from short paths (di- 
rect propagation between the DSPR's involving few or 
no boundary reflections) continue to add coherently for 
large cor, resulting very slow decay. For M| \y most of the 
injected carriers experience a Rashba field nearly aligned 
with their spin. At intermediate Rashba field these yield 



small oscillations that center about a finite value of Rs- 
The other carriers make a contribution to Rs at small 
lor but this becomes incoherent and thus quickly decays 
for large u R (ll|jl|. 

The original idea of the spin transistor involved use of 
an external gate potential, acting in concert with the 
intrinsic confinement potential, to control of the spin 
precession rate We note, however, that gate tun- 
ing of a so for electrons has been experimentally demon- 
strated in relatively few narrow gap semiconductor het- 
erostructures. Two such systems are InGaAs/InP and 
InGaAs/InAlAs In the latter, tuning over about 

a 30% range has been reported. In Fig. || we show how 
this range of tunability translates into a direct modu- 
lation of Rs, for three device widths. Our calculations 
clearly illustrate that the "conventional" spin transistor 
configuration, M||y, (which is most easily fabricated) is 
not optimal — even for a very short channel (L ~ is). 
We find that tunability is maximized for M||x. 

The spin transistor was originally envisaged as a one- 
dimensional device, with only a single populated trans- 
verse subband. Realizable devices in the near term will 
more likely be two-dimensional or, perhaps, quasi-one di- 
mensional channels. Their increased phase space for scat- 
tering can lead to quick suppression of Rs, especially in 
the presence of moderate scattering (TTJ . Hence it ap- 
pears that an extremely narrow channel is a basic re- 
quirement for a spin transistor. 

Our calculations of Rs clarify the important, and 
unique, signatures of spin injected transport in an elec- 
tron gas within a semiconductor. They also point out 
crucial experimental challenges that must be faced in 
making a spin transistor. 
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